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Abstract. In order to clarify the physics of the crossover from a Peierls band insulator to a correlated 
Mott-Hubbard insulator, we analyze ground-state and spectral properties of the one-dimensional half- 
filled Holstein-Hubbard model using quasi-exact numerical techniques. In the adiabatic limit the transition 
is connected to the band to Mott insulator transition of the ionic Hubbard model. Depending on the 
strengths of the electron-phonon coupling and the Hubbard interaction the transition is either first order 
or evolves continuously across an intermediate phase with finite spin, charge, and optical excitation gaps. 

PACS. 71.10.Hf , 71.10.Fd, 71.45.Lr, 71.30.+h 



1 Introduction 

In quasi-one-dimensional materials like halogen-bridged 
: transition metal chain complexes, conjugated polymers, 
organic charge transfer salts, or inorganic blue bronzes 
the itinerancy of the electrons strongly competes with 
electron-electron and electron-phonon interactions, which 
tend to localize the charge carriers by establishing com- 
' mensurate spin- (SDW) or charge-density-wave (CDW) 
\ ground states (GSs). At half-filling, Peierls (PI) or Mott 
(MI) insulating phases are favored over the metallic state. 
Quantum phase transitions between the insulating phases 
are possible and the character of the electronic excita- 
tion spectra reflects the properties of the different insulat- 
ing GSs. A controversial issue is the nature of the PI-MI 
transition and whether or not only one quantum criti- 
cal point separates the PI and MI phases in purely elec- 
tronic model Hamiltonians jl],||,||,[|,|]]. Phonon dynamical 
effects, which are known to be particularly important in 
low-dimensional materials may further modify the 

■ transition. 

In this paper we study the PI-MI quantum phase tran- 
sition in the Holstein-Hubbard model (HHM) at half- 
filling. Exact numerical methods || are used to diagonal- 
ize the HHM on finite chains, preserving the full dynamics 
of the phonons, and the density matrix renormalization 
group (DMRG) technique |)| is applied to the adiabatic 
HHM and the ionic Hubbard model. On finite periodic 
chains we identify one critical PI-MI transition point in 
the HHM where the site-parity of the GS changes and 
the excitation gap in the optical conductivity closes. In 
the adiabatic limit two scenarios emerge with a discon- 
tinuous transition at strong coupling and two subsequent 
continuous transitions in the weak coupling regime with 



the possibility for an intermediate insulating phase with 
finite spin, charge, and optical excitation gaps. 



2 Theoretical models 

The paradigm for correlated electron-phonon systems has 
usually been the one-dimensional HHM defined by 

H = H t -u - g^o J2( b t + +UJ oJ2 b t b i ' C 1 ) 



H,- 



t-u 



-*£( 



cL c 



H.c.) + U^n^nn . (2) 



H t -u constitutes the conventional Hubbard Hamiltonian 
with hopping amplitude t and on-site Coulomb repul- 



sion strength U ; 
nier site i and m. 



creates a spin- a electron at Wan- 
In (mi, the second term cou- 



pies the electrons locally to a phonon created by b\ . Here 
g = ^Je p /ujQ is a dimensionless electron-phonon coupling 
constant, where e p and luq denote the polaron binding en- 
ergy and the frequency of the optical phonon mode, re- 
spectively. 

The GS of the Holstcin model for U = is a Peierls dis- 
torted state with staggered charge order in the adiabatic 
limit ujq — > for any finite e p . As in the Holstein model 
of spinless fermions jic| , p!l| , quantum phonon fluctuations 
destroy the Peierls state for small electron-phonon inter- 
action strength || - an issue which has remained unre- 
solved in early studies of the Holstein model using Monte 
Carlo techniques Above a critical threshold g c (uj ), 
the Holstein model describes a PI with equal spin and 
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charge excitation gaps - the characteristic feature of a 
band insulator (BI). 

The adiabatic limit of the HHM takes the form 

(termed AHHM); it includes the elastic energy of a har- 
monic lattice with a "stiffness constant" K . In this frozen 
phonon approach, A\ = (— l) l A is a measure of the static, 
staggered density modulations of the PI phase. Eq. (||) 
with K — and fixed A is known as the ionic Hubbard 
model (IHM) for which an insulator-insulator transition 
was already established before, although with controver- 
sial results regarding the possibility of an additional in- 
tervening phase Interestingly, the IHM was mo- 
tivated originally in quite different contexts, i.e. for the 
description of the neutral to ionic transition in charge 
transfer salts p3| and ferro-electricity in transition metal 
oxides @,|||. 

3 Numerical results 

3.1 Charge- and spin-structure factors 




0.0 1.0 2.0 



U/2e p 



In order to establish the GS properties of the above models 
and the existence of the PI-MI transition we start with 
the evaluation of the staggered charge- and spin-structure 
factors S c (w) and S s (ir), respectively, 

S ^ = jf £(-1)^^) , Sf = i(n iT - mi) ■ 
j 

Results for the [/-dependences of S c (ir) and S s (n) on an 8- 
site HHM ring are shown in Fig. |l| for two different phonon 
frequencies, corresponding to adiabatic and non-adiabatic 
regimes. The PI regime is characterized by a large (small) 
charge- (spin-) structure factor. Also shown in Fig. [j] 
are results for the AHHM [N = 8 (Lanczos), N = 64 
(DMRG)] which will be discussed in Sec. 3.3 Increasing 
U at fixed e p and uio, Peierls CDW order is suppressed 
as becomes manifest from the rapid drop of S c (tt) which 
decreases nearly linearly in the adiabatic regime, but its 
initial decrease is significantly weaker for higher phonon 
frequencies. The disappearance of the charge ordering sig- 
nal is accompanied by a steep rise in S s (jr) indicating 
enhanced antiferromagnetic correlations in the MI phase. 
The data for S c (ir) provide evidence for a critical point 
U c at which the CDW order disappears: rather abruptly 
for adiabatic and smoothly for non-adiabatic phonon fre- 
quencies. Above U c the low-energy physics of the system 
is qualitatively similar to the pure Hubbard chain; it is 
governed by gapless spin and massive charge excitations. 
As compared to the PI phase, the local magnetic moment 
Li(U/t,e p /t,uJo/t) oc ((Sf) 2 ) is strongly enhanced (e.g. 
L < (8 J 2,l)/L i (0,2,l)~16.4). 



Fig. 1. Staggered charge- (upper) and spin-structure factors 
(middle panel) vs. the rescaled Hubbard interaction U/2e p . 
The lowest panel displays the [/-dependence of the kinetic en- 
ergy Ekin- Lanczos results for the HHM on an 8-site ring are 
given in the adiabatic (triangles) and non-adiabatic (squares) 
regimes. Lanczos (N = 8 ring, down-triangle) and DMRG 
(open 64-site chain, stars) results are shown for the AHHM 
with K — 0.74. Open (closed) symbols belong to GSs with 
site-parity P = —1 (+1). 



We emphasize the weak finite-size dependence of the 
exact diagonalization data for S c (ir) in both the strong- 
CDW Peierls and MI phases (see the upper panel of 
Fig. |^). The variation of S s (n) with the system size N 
points to a vanishing spin-structure factor in the CDW 
state, i.e. for U < U c but, of course, it is beyond our 
current numerical capabilities to perform a real finite-size 
analysis for the HHM with dynamical phonons. In the adi- 
abatic limit, we expect a finite S s (ir) in the MI phase as 
for the corresponding SDW state of the so-called extended 
Hubbard model with U > 2V 

To discuss the lattice dynamical effects in some more 
detail we show in Fig. || the weights of m-phonon states 
in the GS of the HHM at different interaction strengths. 
First of all Fig. || demonstrates that our phonon Hilbert 
space truncation procedure is well-controlled in the sense 
that states with larger numbers of phonons, as accounted 
for in the calculations, have negligible spectral weight. Of 
course, the number of phonons which have to be taken 
into account strongly depends on the physical situation. 
Whereas the GS of the MI is almost a zero-phonon state, 
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of phonons involved in the creation of the Peierls-distorted 
CDW is rather small. On the contrary, for large phonon 
frequencies, g > Q c {oJq) implies A ^> 1 and we observe a 
multi-phonon GS (cf. Fig. y). As a consequence the elec- 
trons are heavily dressed by phonons, forming bipolarons 
in real space, which lower their energy by ordering in a 
staggered CDW. Therefore, in the non-adiabatic strong 
electron-phonon coupling regime, the system is classified 
rather as a charge-ordered bipolaronic insulator than as a 
BI. The kinetic energy is much more suppressed for the 
bipolaronic than for the band PI (cf. lower panel of Fig. [l]). 
Since SDW correlations reduce Ek in as well, the kinetic 
energy reaches a maximum when the system crosses from 
the PI to the MI regime. 
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Fig. 2. Finite-size scaling of charge- (upper panel) and spin- 
(lower panel) structure factors in the HHM. 



multi-phonon states become increasingly important if U 
is reduced (e p is enhanced) in the PI regime. 

For small U and low phonon frequencies the PI phase 
appears for g > g c (u>o), even if the ratio A = s p /2t is small. 
In the GS of such a conventional BI phonons populate pre- 
dominantly the q — and n modes, but the total number 
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Fig. 3. Phonon distribution in the GS of the HHM for various 
model parameters. In the MI state (open symbols) the weight 
of the zero-phonon state is almost one, |c m | 2 ~ 1. 



Valuable insight into the nature of the PI-MI transition is 
obtained from symmetry considerations The BI-MI 

transition of the IHM on finite lattices was shown to be 
connected to a GS level crossing with a site-parity change, 
where the site inversion symmetry operator P is defined 
by Pc4.Pt = c ] v -_ i(7 with N = An for i = 0, ...,7V - 1. 
This feature will become evident in the regular part of 
the optical conductivity at T = 0, 



a reg ( ) = £. I (V'O |5' I V'm) | 2 

[ ' N 2-? E, 



E 



S{uj-E m + E Q ). (4) 



Here, IV'o) and \ip m ) denote the GS and excited states, re- 
spectively, and E m the corresponding eigenenergies. Im- 
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Fig. 4. Optical conductivity in the 8-site HHM for ujo = 0.1t 
and g 2 = 7- Top panel: PI phase for U = 0; middle panel: near 
criticality U ~ U pt', lower panel: MI phase for U = 3t. Dashed 
lines give the normalized integrated spectral weights S re9 (u). 
The lower two panels include a res for g — (dotted lines), i.e. 
for the pure Hubbard chain. 
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portantly, the current operator j — — iet^io-( c L c i+icr — 
c| +1 a c ia ) has finite matrix elements between states of dif- 
ferent site-parity only. 

The evolution of the frequency dependence of a reg {uo) 
from the PI to the MI phase with increasing U is illus- 
trated in Fig. [|. In the PI regime the electronic exci- 
tations are gapped due to the pronounced CDW corre- 
lations. The broad optical absorption band for U = 
results from particle-hole excitations across the BI gap 
which are accompanied by multi-phonon absorption and 
emission processes. The shape of the absorption band re- 
flects the phonon distribution function in the GS. Exci- 
tonic gap states may occur in the process of structural 
relaxation. At U op t the optical gap A opt closes, and due 
to the selection rules for optical transitions this necessar- 
ily implies a GS level crossing with a site-parity change. 
We have explicitly verified that the GS site parity in the 
PI phase is P = +1 and P = — 1 in the MI phase (see also 
Fig. |]). For the HHM on finite rings U op t is identical to 
the critical point where S c (tt) sharply drops. 

For the adiabatic phonon frequency used in Fig. ^| the 
phonon absorption threshold is small and, because the GS 
is a multi-phonon state, we find a gradual linear rise of 
the integrated spectral weight S reg (uj) = f Q a re9 (uu') du' . 
S reg (oj) j S reg (oo) is a natural measure for the relative 
weight of the different optical absorption processes. In 
contrast, in the non- adiabatic regime (loq > t), the lowest 
optical excitations have mainly pure electronic character 
in the vicinity of U op t- As a result the gap is closed by a 
state having large electronic spectral weight. 

In the MI phase the optical gap is by its nature a cor- 
relation gap. The lower panel in Fig. [| shows clearly that 
a(cu) of the HHM in the MI phase is dominated by excita- 
tions which can be related to those of the pure Hubbard 
model. In addition, phononic sidebands with low spectral 
weight and phonon-induced gap states appear. 



3.3 Phase diagram in the adiabatic limit 

The above results for the HHM establish the PI-MI phase 
transition scenario on small rings and trace it to the level 
crossing of the two site-parity sectors. In order to draw 
conclusions about the phase diagram in the adiabatic 
regime we exploit the connection to the AHHM. 

The magnitude of S c {n) in the HHM for [7 = 
and luq = O.li allows a straightforward way to fix the 
stiffness constant K in Eq. (g). Using the result of the 
AHHM for S c {ir) at U = K = we determine first 
the ionic potential strength Aq by the requirement that 
S I c HM {ir,A ) = S? HM (%) for the same chain length and 
periodic boundary conditions. In a second step, the GS en- 
ergy of the AHHM, E (K, A, U = 0), determines K by the 
criterion that Eo is minimized for A = Aq. We thereby 
obtain K = 0.74, which is henceforth kept fixed when 
the interaction U is turned on. For each value of U, the 
ionic potential strength of the AHHM is then obtained 
by minimizing E$(K, A,U) with respect to A, yielding 
A = A(U, K) as shown in Fig. || (triangles). The resulting 
structure factors for the AHHM are plotted in Fig. |l|, too, 




U/2e 

p 

Fig. 5. Level crossing line A cr (U) of the IHM for an 8-site ring 
(diamonds) and from extrapolating Lanczos data for N < 14 
to a 64-site chain (circles). In addition: ionic potential strength 
A(U, K) of the AHHM for an 8-site ring (triangles) and on an 
open 64-site chain (stars, DMRG results) for K — 0.74. 



and agree very accurately with the 8-site HHM ring data 
for uo/t = 0.1. This agreement reconfirms numerically 
that the AHHM is indeed the appropriate effective model 
to describe the CDW phase of the HHM in the adiabatic 
limit. The drop in S c (ir) at the transition point results 
from a discontinuous vanishing of A(U, K) (see Fig. ||). 
The large charge structure factor S c (tt) below U c and the 
enhancement of the spin structure factor S s (n) above U c 
as well as the sharp changes at the transition point find a 
natural explanation with the results for A(U, K) in Fig. [|. 
Below the transition A is finite implying long range CDW 
order in the GS. At the transition point A vanishes dis- 
continuously and thereby the AHHM reduces to the pure 
Hubbard model (A = 0). 

Given the value for the stiffness constant K we also 
plot in Fig. | A(U, K) obtained fr om DMRG on an open 
chain of length N = 64. For comparison, the correspond- 
ing results for S c (tt) and S s (ir) in the AHHM are shown in 
Fig. [l, too (stars). S c (ir) decreases smoothly and almost 
linearly; although unresolved on the vertical scale in Fig. [j] 
the transition remains discontinuous as a consequence of 
the results for A(U, K) in Fig. M. In contrast to the be- 
havior of the 8-site chain, A(U, K) here decreases more 
smoothly with increasing U and vanishes discontinuously 
near U/2e p « 0.75. The small discontinuous increase in 
S s (ir) at the transition is also hardly resolved for the 64- 
site chain in contrast to the 8-site chain data. The discon- 
tinuous nature of the PI-MI transition in the AHHM for 
Ai = (— l) l A is obvious in the atomic limit t = where 
A = 1/K for U < U c = 1/K and A = for U > U c . 
As verified above, the first order nature persists for finite 
small t, i.e. in the strong coupling regime U, K^ 1 3> t. 

Also shown in Fig. H is the level crossing line A cr (U) 
of the IHM for N = 8 (diamonds) and N = 64 (circles) 
chain. A cr (U) for N = 64 was obtained from extrapolat- 
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Fig. 6. Insets (1) - (5): Evolution of the ground-state energy 
vs. A in the AHHM in different regions of the (K ,U) pa- 
rameter plane. From the variations in E(A) a crossover from 
a discontinuous PI (with A > 0) to MI (A = 0) transition 
to a second order transition is deduced. Main figure: Phase 
diagram of the AHHM; the solid line represents a discontinu- 
ous, first order and the dashed line a continuous second order 
transition. These results summarize Lanczos data for a 14-site 
AHHM chain with periodic or open boundary conditions. De- 
tailed runs where performed for U = 0.3t and U = 5t. A possi- 
ble additional continuous transition (dotted line) between two 
insulating phases with finite A is indicated as well. 



ing Lanczos results for rings of up to 14 sites to a 64-site 
chain jl8| . Importantly, A(U, K) and A cr (U) do not inter- 
cept because A(U,K) jumps to zero before reaching the 
level crossing point of the IHM. 

The DMRG results presented in Fig. | for N = 64 raise 
the question whether the discontinuous transition in the 
AHHM can turn into a continuous transition on approach- 
ing the weak coupling regime by increasing the stiffness 
constant K. Indeed, as we have explicitely verified by ex- 
act diagonalization of a periodic (and open too) AHHM 
ring of length N = 14, the transition is second order in the 
regime [/, K^ 1 <C t. The corresponding Lanczos results 
for the variation of the GS energy vs. A in the (K, U)- 
parameter plane are summarized in Fig. ©. Detailed K- 
scans were performed for weak (U — 0.3i) and strong 
(U = 5t) Hubbard interaction. The evolution of E(A) 
in the AHHM in fact reveals that the transition from the 
PI to the MI phase (sequence (2) — (3) — (4)) occurs dis- 
continuously at strong coupling K~ 1 ,U ^> t, while the 
transition follows a Ginzburg-Landau-type behavior for a 
second order phase transition at weak coupling (sequence 
(1)- (5) in Fig.©). 

Due to the continuous decrease of A(U, K) at weak 
coupling AjU, K) necessarily intercepts the A cr (U) line of 
the IHM |l8|. This intercept marks the point U op t when 
the site-parity sectors become degenerate and the optical 
absorption gap A opt disappears. This situation therefore 
implies the existence of an intermediate region U op t < U op t < U < U s in which all excitation gaps arc finite. The 
U < U s with finite A, where U s marks the point where A CDW persists for all U < U s . The site-parity eigenvalue 
continuously vanishes. Since A = for U > Us, i.e. when is P = +1 in the PI and P = — 1 in the MI phase. 



Fig. 7. Qualitative behavior of the excitation gaps versus U 
in the weak coupling regime of the AHHM U, K' 1 < t. Solid 
line: optical excitation gap A opt , dotted line: charge gap A c , 
dashed line: spin gap A 3 . PI phase: A c — A s and site parity 
P = +1; MI phase: A opt = A c and P = -1. 



the AHHM reduces to the Hubbard model, the spin gap 
vanishes at U s . The intermediate insulating phase thus has 
finite spin, charge, and optical excitation gaps. For weak 
coupling the PI-MI transition therefore evolves across two 
critical points U op t and U s . The U vs. K" 1 phase diagram 
contains a multicritical point at which a first order line 
splits into two continuous transition lines. The additional 
transition line is also indicated in Fig. © (dotted line). 
For weak U at fixed K the transition at U op t is expected 
to be of Kosterlitz-Thouless type since it corresponds to 
the merging of the energies of the two site-parity sectors; 
the CDW vanishes in a second order type transition at 
U = U s . The E(A) behavior, however, can only detect 
the boundary to the MI phase of the AHHM where the 
GS energy is minimized for vanishing A. 

We summarize these findings in the diagram for the 
excitation gaps shown in Fig. [?], 



A c = E (N/2 + l,N/2) + E (N/2 - 1, N/2) 

-2E (N/2,N/2), 
A s = E a (N/2 + 1, N/2 - 1) - E (N/2, N/2) 



(5) 
(6) 

where Eq (iV-f , iVjJ is the GS energy of the system with 
spin- up and N± spin-down electrons. In the Peierls 
BI phase for U < U op t the spin and charge gaps, are 
equal and finite and remarkably A opt ^ A c (for a sim- 
ilar conclusion in the IHM see |Q). At U = U op t the 
site-parity sectors become degenerate, A opt — but re- 
markably A c — A s > 0. For U > U s the usual MI phase 
of the half-filled Hubbard chain with A opt = A c > A s = 
is realized. For strong coupling, when the PI to MI transi- 
tion is first order, U opt = U s , the spin gap discontinuously 
disappears at the transition and the optical gap jumps 
from zero to the finite charge gap value of the Hubbard 
chain. In weak coupling there exists an intermediate region 



6 



H. Fehske et al.: Nature of the Peierls- to Mott-insulator transition in ID 



The insulating, intermediate phase at weak coupling as 
identified above remains yet to be characterized. For the 
insulator-insulator phase transition(s) in the IHM Fab- 
rizio et al. proposed the existence of an intermediate phase 
with a long range bond order wave (BOW) based on a 
bosonization analysis B. BOW order is characterized by 
a finite expectation value of the staggered bond charge 

B = lEwf-il'^/.+i^ + Hx -)- Some Positive nu- 
merical evidence has indeed been reported for enhanced 
BOW correlations above the level crossing transition in 
the IHM Yet, these results have remained am- 

biguous so far and no consensus has been reached about 
the existence of long range BOW order in the IHM. In 
the attempt to search for BOW order the DMRG calcula- 
tions, which for numerical accuracy reasons are predom- 
inantly performed on open chains, suffer from the fact 
that Friedel-like bond charge density oscillations are in- 
duced by the chain ends already for the pure Hubbard 
chain (|). The identification of BOW order in the IHM 
or AHHM by DMRG on open chains therefore requires a 
delicate subtraction procedure to discriminate a BOW sig- 
nal from the edge induced bond charge oscillations of the 
Hubbard chain. We have nevertheless attempted to search 
for BOW correlations in the weak coupling regime of the 
AHHM, where the continuous nature of the transition into 
the MI phase was established by the Lanczos results on 
the periodic or open 14-sites chain, i.e. these calculations 
naturally focused on the weak-[/ regime (U < t). This 
weak coupling regime is notoriously hard for numerical 
evaluations; unfortunately the numerical accuracy needed 
to allow a firm conclusion about the presence or absence 
of a BOW signal could not be achieved within our DMRG 
runs. 

While a confirmation is thus still lacking BOW or- 
der remains a vivid candidate order in coexistence with 
a CDW to characterize the intermediate phase in the 
AHHM at weak coupling. We furthermore note that if the 
existence of a BOW is verified in the AHHM, its phase di- 
agram would be remarkably similar to the extended Hub- 
bard model with nearest neighbor Coulomb repulsion with 
an intervening BOW phase in the crossover between the 
CDW and MI phases at weak coupling p2|,p3[. 



4 Conclusions 



In summary, we have found a PI-MI transition in the HHM 
above a threshold electron-phonon coupling. The transi- 
tion results from a GS level crossing with a change in the 
GS site-parity eigenvalue. In the adiabatic limit two sce- 
narios emerge with a discontinuous PI-MI transition for 
U, K^ 1 3> i, and two continuous transitions for weak cou- 
pling U, K^ 1 <C t with an intermediate phase where CDW 
order persists. In the non-adiabatic regime our structure 
factor data indicate that the PI-MI transition proceeds 
continuously 
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